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ABSTRACT 

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H- ATLAS SDP observations, to 
investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these 
models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of 
two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good 
sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we 
support highly obscured star formation that contributes little to the optical extinction as a more probable explanation. 
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1. Introduction 

Edge-on spiral galaxies are an important class of galaxies in 
which the distribution and properties of interstellar dust grains 
can be studied in great detail. The dust in these systems shows 
prominently as dust lanes in optical images; for several edge- 
on spiral galaxies, the dust distribution has been modelled by 
fitting realistic radiative transfer models to such optical images 
(Kylafis & Bahcall 1987; Xilouris et al. 1997, 1998, 1999; Alton 
et al. 2004; Bianchi 2007). The conclusion of these works is that, 
in general, the dust disc is thinner (vertically) but radially more 
extended than the stellar disc and that the central optical depth 
perpendicular to the disc is less than one in optical wavebands, 
making the disc almost transparent when seen face-on. 

A complementary way of studying the dust component in 
galaxies is through its thermal emission at far-infrared (FIR) 
and submm wavelengths. A self-consistent treatment of extinc- 
tion and thermal emission, i.e. a study of the dust energy bal- 
ance, gives the strongest constraints on the dust content of spi- 
ral galaxies. When we quantitatively compare the results from 
extinction studies and FIR/submm emission studies, a discrep- 
ancy is found: the relatively optically thin dust discs emerging 
from the optical modelling absorb about 10% of the stellar ra- 
diation, whereas FIR studies of normal spiral galaxies indicate 
they typically reprocess about 30% of the UV/optical radiation 
(Fopescu & Tuffs 2002). This problem is clearly illustrated when 
applied to individual edge-on spiral galaxies: the predicted FIR 
fluxes of self-consistent radiative transfer models that success- 
fully explain the optical extinction generally underestimate the 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led principal- investigator consortia and with impor- 
tant participation from NASA. 



observed FIR fluxes by a factor of about three (Fopescu et al. 
2000; Misiriotis et al. 2001; Alton et al. 2004; Dasyra et al. 
2005). Several scenarios (see Sect. 4) have been proposed to ex- 
plain this discrepancy, but a major problem distinguishing be- 
tween these is that the number of edge-on galaxies for which 
such detailed studies have been done so far is limited, owing to 
the poor sensitivity, spatial resolution, and limited wavelength 
coverage of the available FIR instruments. Herschel (Pilbratt 
et al. 2010) offers the possibility of studying the dust energy- 
balance in spiral galaxies in more detail than before. The com- 
bination of the sensitivity and the wavelength coverage of the 
PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) in- 
struments, together covering the 70 to 500 jim region where the 
emission from cold dust dominates, enables us to make reliable 
estimates of the total thermal emission of the interstellar dust. 

A modest number of large, nearby edge-on spiral galax- 
ies are targetted by Herschel key programs, including the Very 
Nearby Galaxies program and the Herschel Reference Survey 
(Boselli et al. 2010). Much promise, however, is offered by 
the Herschel Astrophysical TeraHertz Large Area Survey (H- 
ATLAS, Eales et al. 2010), which will image a total area 
of 550 deg2 with PACS and SPIRE at 100, 160, 250, 350, 
and 500 /xm. Based on the Revised Flat Galaxy Catalogue 
(Karachentsev et al. 1999) and the 2MASS-selected Flat Galaxy 
Catalogue (Mitronova et al. 2004), we expect between 100 and 
200 resolved edge-on spiral galaxies in the H-ATLAS survey 
area. 

In this Letter, we investigate the dust energy balance in the 
nearby spiral galaxy UGC 4754 as a first case. This Sbc galaxy is 
the nearest (v,. = 5662 ± 45 km/s; Jones et al. 2009) and largest 
edge-on spiral galaxy in the 16 deg^ field observed during the 
H-ATLAS science demonstration phase (SDP). We present de- 
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Fig. 1. The multi- wavelength view of UGC 4754. The top row shows the Herschel images, the bottom row a selection of optical and 
NIR images from SDSS and UKIDSS. The field-of-view of each image is 100" x 100". 



tailed radiative transfer simulations of this galaxy that fit optical 
and near-infrared imaging data, compute synthetic FIR fluxes 
and images, and compare the model predictions to the Herschel 
observations at 100, 160, 250, 350, and 500 yum. In Sect. 2 we 
present the observations and data reduction, in Sect. 3 we present 
our radiative transfer modelling and compare the model results 
to the observations, and in Sect. 4 we discuss these results and 
present out conclusions. 

2. Observations and data reduction 

The edge-on spiral galaxy UGC 4754 was observed on 22 
November 2009 with the PACS and SPIRE instraments on- 
board Herschel as part of the H-ATLAS SDP obsei-vations of 
a 4 X 4 deg2 equatorial field centred on (o-, S) ~ (9^05™, +0°30'). 
Data were gathered simultaneously in the green and red PACS 
bands (100 and 160 //m) and the three SPIRE bands (250, 350, 
and 500 pm\). The PACS and SPIRE time-line data were reduced 
using HIPE, with reduction scripts based on the standard reduc- 
tion pipelines. Maps from the SPIRE data were produced using 
a naive mapping technique after removing the eff'ects of temper- 
ature variation on the timelines using a method developed by the 
H-ATLAS consortium (Pascale et al. 2010). Noise maps were 
generated by using the two cross-scan measurements to estimate 
the noise per detector pass, and then for each pixel the noise is 
scaled by the square-root of the number of detector passes. Maps 
from the PACS data were produced using photproject (Ibar et al. 
2010). 

UGC 4754 was clearly detected and spatially resolved in all 
the observed bands (Rigby et al. 2010, in prep.; Smith et al. 
2010a, in prep.). The total fluxes are 894 ± 182 mJy at 100 /im, 
1287±209 mJy at 160 yum, 878+1 34 mJy at250jum, 383+61 mJy 
at 350 jum, and 122 ±20 mJy at 500 //m. The quoted errors in- 
clude statistical errors, confusion noise, and a 15% absolute cal- 
ibration uncertainty (Poglitsch et al. 2010; Griffin et al. 2010). 

Apart from the Herschel images, we used optical ugriz im- 
ages from SDSS DR7 (Abazajian et al. 2009) and near-infrared 
YJHK images from UKIDSS DR5 (Lawrence et al. 2007) for 
our analysis. These images were obtained from the SDSS and 
UKIDSS archives and reduced and flux-calibrated using stan- 
dard reduction methods in MIDAS. Total fluxes were measured 



from integrating over the flux-calibrated maps after removal of 
the foreground stars. Finally, we also used IRAS fluxes at 60 and 
100 /im from the Faint Source Catalog (Moshir et al. 1990). 

3. Analysis 

The first step in our analysis was to construct a radiative trans- 
fer model for UGC 4574, based on the gri and YJHK images; 
the SDSS u and z bands were not used in our fitting procedure 
because of low signal-to-noise. We used the 3D Monte Carlo ra- 
diative transfer code SKIRT for the modelling. This code was 
initially developed to investigate the effects of dust extinction on 
the photometry and kinematics of galaxies (Baes & Dejonghe 
2002; Baes et al. 2003), but has evolved to a flexible radiative 
transfer tool that can model the absorption, scattering, and ther- 
mal emission of circumstellar discs and dusty galaxies (e.g. Baes 
et al. 2005; Gadotti et al. 2010). 

The stellar distribution was represented as a combination 
of a double exponential disc (i.e. a disc which follows an ex- 
ponential profile both radially and vertically) and a flattened 
Sersic model. The intrinsic stellar SED of the galaxy' was de- 
termined by fitting a stellar population synthesis model to the 
dereddened SDSS-UKIDSS flux densities and adapting the code 
used in Hatziminaoglou et al. (2008, 2009); the resulting model 
SED corresponds to a population of 8 Gyr old with an expo- 
nentially decaying star formation rate and an initial burst dura- 
tion of 0.15 Gyr. The dust component was also represented as 
a double exponential disc. The optical properties of the dust are 
based on the BARE_GR_S model of Zubko et al. (2004), which 
accurately reproduces the extinction, emission, and abundances 
in the MiUcy Way. Our final model contains 11 free parameters: 
the stellar scalelength, stellar scaleheight, bulge effective radius, 
Sersic index, bulge flattening, bulge-to-disc ratio, total luminos- 
ity, dust scalelength, dust scaleheight, total dust mass (or equiv- 
alently, optical depth), and inclination. 



' To limit the numbers of parameters, we assumed a single intrinsic 
SED for the entire galaxy, i.e. we assumed no intrinsic colour gradients 
and the same SED for the bulge and disc. This assumption has virtually 
no effect on the predicted FIR emission, which is the main goal of this 
Letter. 
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Fig. 2. Comparison of the observed r, J and K band images 
of UGC 4754 (top row) and the resuhs of our SKIRT radiative 
transfer modelHng (bottom row). All images show only the cen- 
tral 30" x 30". 



For every choice of these 1 1 parameters, a set of synthetic 
images can be constructed by solving the radiative transfer equa- 
tion (taking both absorption and scattering into account), con- 
volving the resulting radiation field with the appropriate trans- 
mission curves and PSF. Using a x^ minimization, we looked 
for the set of parameters for which the coiTesponding synthetic 
images reproduced the obsei-ved images best. The initial input 
values for the parameters of the stellar distribution were deter- 
mined from a bulge-disc decomposition of the /-band image us- 
ing the BUDDA code (de Souza et al. 2004; Gadotti 2008). For 
the dust distribution, the initial input values were calculated us- 
ing the average dust-star geometry from Xilouris et al. (1999). 
In general, our fitting approach is very similar to the method 
applied by Xilouris et al. (1999) and Bianchi (2007). The main 
difference is that we fit a single photometric model to all optical 
and NIR images simultaneously, whereas Xilouris et al. (1999) 
and Bianchi (2007) fit a new model to each different band. 

Figure 2 shows the final radiative transfer model compared 
to the observed images in the r, J, and /T bands. The stellar disc 
in this model has a scalelength of 4.05 kpc and a scaleheight of 
330 pc, whereas the dust has a scalelength of 6.1 kpc and a scale- 
height of 270 pc. The total bolometric luminosity of the galaxy 
is 1.8 X 10'*' Lq, of which the bulge only contributes 8%. The 
total dust mass in our best- fitting model is 1 .0 x 10^ M^ and this 
translates into face-on optical depths (measured along the entire 
z-axis) of 0.73, 0.49, 0.17 and 0.07 in the g, r, 7, and /Tbands, re- 
spectively. The corresponding edge-on optical depths (measured 
along the entire line-of-sight through the centre) are 16.5, 11.1, 
3.8, and 1.6. In our model, 2.1 x lO'* Lq or 12% of the bolo- 
metric stellar radiation is absorbed by the dust. Both the relative 
star-dust geometry, with a naiTow and extended dust distribution, 
and the modest values of the optical depth agree with similar ra- 
diative transfer modelling results of other edge-on spiral galaxies 
(Xilouris et al. 1999; Alton et al. 2004; Bianchi 2007). 

The next step in our analysis is to compare the resulting 
FIR emission of our radiative transfer model to the Herschel 
observations. Calculation of the FIR emission of a radiative 
transfer model is carried out using SKIRT. At each position in 
the galaxy, the mean intensity of the radiation field is calcu- 
lated at every wavelength during the radiative transfer simula- 
tion. From this mean intensity, the equilibrium dust tempera- 
ture of each species of dust grains can be calculated using the 
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Fig. 3. The optical/NIR (top) and FIR/submm (bottom) spec- 
tral energy distribution of UGC 4754. See label and text for the 
meaning of the different lines. 

energy balance equation. We considered 10 different size bins 
for each of the three dust materials in the Zubko et al. (2004) 
model (graphite, silicate, and PAHs). We did not consider tran- 
sient heating by very small grains and PAHs in our simulation, 
but to estimate the impact of transiently heated grains on the 
SED in the Herschel region, we calculated the FIR emission in 
two ways. First, we assumed that all grains are in LTE and we 
added all of their thermal emission to calculate the FIR emis- 
sion. Consequently, we repeated this, but only took the emission 
of big grains {a > 0.01 ;um) into account. 

Figure 3 shows the spectral energy distribution (SED) of 
UGC 4754 at optical/NIR (top panel) and FIR/submm (bottom 
panel) wavelengths. The solid black line in the top panel corre- 
sponds to the attenuated SED of the SKIRT model fitted to the 
SDSS and UKIDSS images, and the magenta line is the unat- 
tenuated SED. Our radiative transfer model reproduces the opti- 
cal/NIR SED very well. In the bottom panel, the solid black line 
coiTesponds to the FIR emission of the model assuming LTE for 
all grains, and the dotted line represents the contribution of only 
the large grains. Our model significantly underestimates the ob- 
served IRAS and Herschel fluxes. This behaviour, a dust energy- 
balance problem, has been noted for other spiral galaxies where 
the FIR emission calculated from radiative transfer modelling 
is compared to the observed IRAS, ISO, and/or SCUBA fluxes 
(e.g. Bianchi et al. 2000; Popescu et al. 2000; Misiriotis et al. 
2001). 

4. Discussion and conclusion 

An apparently straightforward explanation to the dust energy- 
balance problem might seem that the standard FIR emissivity, 
derived from dust models finetuned to the diffuse dust emission 
in the Milky Way, is actually underestimated by a factor two 
to three. This idea was advocated by Alton et al. (2000, 2004) 
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and Dasyra et al. (2005) to explain the excess at SCUBA wave- 
lengths of several edge-on spiral galaxies. Support for this ex- 
planation is the wide range of empirical values for the FIR emis- 
sivity that have appeared in the literature (see e.g. Hughes et al. 
1997; Alton et al. 2004). Thanks to Herschel's wavelength cov- 
erage between 100 and 500 yum, we can rule out this possibility, 
because we can accurately compare the absorbed stellar lumi- 
nosity in our radiative transfer model (Labs = 0.12 Lboi) to the 
observed luminosity L^ emitted by the dust. We calculated the 
observed dust luminosity LpiR in the FIR (between 80 //m and 
1 mm) by integrating the SED of the best-fitting modified black- 
body model to the Herschel data. This model, indicated as the 
green line in the bottom panel of Fig. 3, corresponds to a modi- 
fied blackbody model with M^ = 1.86 x 10^ Mq and T = 19.6 K 
and yields LpiR = 5.0 x 10"^ L© = 0.27 Lboi- Note that LpiR is a 
lower limit to L^, since the mid-infrared luminosity (dominated 
by non-equilibrium emission by very small grains and PAHS) 
has not yet been taken into account. These numbers clearly indi- 
cate that there is a real discrepancy in the energy balance. This 
discrepancy cannot be lifted by changing the FIR emissivity at 
some wavelength; indeed, the energy balance equation impies 
that the ratio Labs/^d is independent of the normalization of the 
emissivity. 

The most likely explanation for this energy balance prob- 
lem is that a sizable fraction of the FIR/submm emission arises 
from additional dust that has a negligible extinction on the bulk 
of the starlight. One option is the presence of a second, thinner 
dust disc, originally proposed by Popescu et al. (2000). Such a 
configuration has been successful in explaining the energy bal- 
ance of spiral galaxies (Misiriotis et al. 2001, Popescu et al. 
2010, in prep.) and the observed attenuation-inclination relation 
(Driver et al. 2007), but its validity has been questioned based 
on deep/: band images of NGC 891 (Dasyra etal. 2005). An al- 
ternative hidden dust source are young stars deeply embedded in 
dusty molecular clouds. The compact dust clumps can boost the 
FIR/submm emission of the dust, while keeping the extinction 
relatively unaltered (e.g. Silva et al. 1998; Bianchi et al. 2000; 
Popescu & Tuff^s 2005; Bianchi 2008; Michalowski et al. 2008). 

A first indication that embedded star-forming clouds might 
be the solution to the case of UGC 4574 is that the discrepancy 
between our radiative transfer model and the observed FIR SED 
is stronger at shorter than at longer wavelengths. This implies 
that warmer dust (such as in star-forming regions) is necessary to 
bring the model in balance with the data. Further evidence is the 
forthcoming study by Smith et al. (2010b, in prep.) who analyse 
the SED of 1260 galaxies from the H- ATLAS SDP field using 
the technique of da Cunha et al. (2008). This technique consists 
of fitting the observed SED with a combination of attenuated 
stellar emission, dust emission from star forming-regions, and 
dust emission from the general ISM. For UGC 4754, Smith et al. 
(2010b) find that more than 50% of the bolometric luminosity of 
the galaxy is reprocessed by dust. In particular, their modelling 
predicts a significantly bluer intrinsic UV continuum. The very 
efficient absorption of this radiation in star-forming regions can 
power the FIR luminosity, while almost remaining unnoticed in 
large-scale optical extinction maps. Combining the results from 
that study with ours, we conclude that embedded star formation 
is the most likely way to lift the energy balance in UGC 4754. 

The present case study of UGC 4754 is a demonstration of 
what we can learn from detailed dust energy balance studies of 
edge-on spiral galaxies. In future contributions we will apply 
similar modelling to larger sets of edge-on galaxies and investi- 
gate possible systematic links with other galaxy parameters; for 
example, one would expect a correlation between the strength 



of the dust energy balance problem and the star formation rate 
if embedded star formation is the main driver (e.g. da Cunha et 
al. 2008, 2010). The Herschel observatory and, particularly, the 
H- ATLAS project offer the possibility to perform such studies in 

a systematic way. 
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